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Abstract 

We introduce an eikonal Monte Carlo model (running in con- 
junction with HERWIG) for simulating multiparticle production in 
hadron-hadron interactions. We compare our simulated data to the 
CDF Tevatron measurement of the underlying event activity in hard 
inelastic proton-antiproton scattering at - v /s=1.8 TeV. By fixing the 
only free parameter in our model, the total hadron-hadron cross sec- 
tion, we find that our model describes the data better than either the 
HERWIG Underlying Event model or the Hard Multiparton model. 



1 Introduction 



In this paper we introduce an eikonal Monte Carlo model which simulates 
hadron-hadron collisions in which a "hard" parton-parton scattering has oc- 
curred. In hard events, described by perturbative QCD, the outgoing parti- 
cles with large transverse momentum form jets. These are clusters of outgoing 
hadrons produced approximately in the direction of the two outgoing partons 
of the hard parton-parton scattering. Jet physics plays an essential role as a 
test of perturbative QCD and in determining the parameters of the Standard 
Model. It will be important for future detection of new particles and might 
help shed light on the nature the Pomeron. 

If we are to find some new physics at future (or current colliders) or 
test the known physics (such as perturbative QCD), we will need to give an 
accurate description of the already known physical process which would then 
be considered as a background to the new physics. 

The underlying event, in hard interactions, is defined to be all the addi- 
tional interactions that are not part of the hard scattering. Thus the under- 
lying event may contain particles from proton-antiproton, remnant-remnant 
interactions and hard multiparton interactions, if such interactions do occur. 

At present, our understanding of underlying event physics is not sufficient 
to match the precision of experimental measurements, so to fully exploit 
current and future data we need to improve that understanding. 

To this end the CDF collaboration recently measured the activity of the 
underlying event in proton-antiproton hard scattering [Q]. They found that 
its activity is considerably larger than in soft collisions and that none of 
the QCD Monte-Carlo models describe correctly all the properties of the 
underlying event (however the best results were obtained with Pythia 0). 

Our goal is to introduce a simple eikonal Monte Carlo model, with a 
small number of parameters, and compare it to the measured data from the 
CDF collaboration 0. In Sec. 2 and Sec. 3 we present the overall theoretical 
framework of the model and give a description of the hard part, in Sec. 4 we 
present our assumptions behind the soft subprocess and its implementation 
within the HERWIG Monte Carlo package and in Sec. 5 we compare our 
theoretical results with those from the CDF experiment PJ and two other 
underlying event models, the HERWIG underlying event model || and the 
Hard Multiparton model 
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2 Optical theorem and impact parameter rep- 
resentation 



In this section we will review some well known expressions for cross sections 
within the eikonal approximation, which we will use as the basic framework 
for our model. 

For two incoming particles A, B and two outgoing particles C, D we write 
their four-momenta as Pa,Pb,Pc,Pd and define the invariant Mandelstam 
variables s and t as: 

s = (p A +p B ) 2 , t= (p A -p c ) 2 . (1) 

The optical theorem relates the imaginary part of the elastic scattering 
amplitude A(s, t = 0) to the total cross section a tot through a unitarity 
relation ||: 

a tot = -Im(A(s,t = 0)). (2) 

s 

Neglecting effects that depend on the spin of the scattering particles we can 
express the scattering amplitude as the Fourier transform of the elastic scat- 
tering amplitude in impact parameter a(b, s) (b being the impact parameter 
and q being the transverse component of momentum transfer in the CMS) 
as : 

A(s, t) = 4vrs J db 2 a{b, s)e iq - b . (3) 
Using the optical theorem (2) we can express the cross sections as: 

a tot = An J db 2 Im(a(h, s)), (4) 
(Tela = 4tt y db 2 \a(b, s)\ 2 , (5) 

<?inel = °tot — Vela- (6) 

We express the elastic amplitude in terms of the eikonal function, s), as 
& 

P -X(b,s) _ i 

a(b,s) = . (7) 

With the eikonal expression (7) we can write the cross sections as : 

roc 

a tot = 2n db 2 [l - e-* b >% (8) 
Jo 
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&ela — 



7T / db 2 [1 - e - x(b ' s) ] , 



Jo 



(9) 



/ 

Jo 



oo 



Cinel = 



dfe 2 [l- e - 2 x(M]. 



(10) 



We choose the total cross section in (8) as the main parameter of our model. 
In the next section we will specify the expression for \ip, s). 

3 Expression for the eikonal 

In the following section, we will give an expression for the eikonal xipi s ) m 
(8) and provide some justification for its expression. 

3.1 The hard part of the eikonal 

We start by assuming that at high center of mass energy (y/s > 100 GeV), 
the main hard scattering is accompanied by semihard interactions, which 
arise as a result of hard multiple parton scattering in which each parton 
carries a very small fraction of its parent hadron's momentum. Multiparton 
interactions lead to the appearance of minijets, i.e. jets with transverse 
energy much smaller than the total energy (y/s) available in the hadron- 
hadron collision. Multiparton interactions will then drive the rise of the 
inclusive jet cross section at high energies depending on the small- a; behaviour 
of the parton distribution or, to a good approximation, that of the dominant 
gluon distribution. 

In refs.[0 |§, 0, [|J, 0] it was shown that the average number of sec- 
ondary hard scatters in an event can be deduced from the QCD improved 
parton model. For our hard multiple-interactions we follow the model pre- 
sented in jl where the average number of secondary hard scatters (n(b, s P p)) 
is given by : 



where the profile function A(b) specifies the overlap of partons in the two 
hadrons in impact parameter. It can be written as a convolution of form 
factor distributions of two incoming hadrons, 



(11) 




(12) 



for which we take 



G p (b) = Gp(b) = J 



d 2 h exp (k ■ b) 

(2^)(l + k 2 /^ 2 ) : 



(13) 
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with /i 2 = 0.71 GeV 2 . The integral then yields [§: 

A(b) = j£-b*)*KsQjb), (14) 
yo7r 

with Ki(x) the modified Bessel function. The overlap function A(b) satisfies 

1. (15) 



J 7idb 2 A(b) 



The cross section appearing in (11), cr# c (s P p), is the inclusive cross section 
for pp^ jets with p t >Ptmin ■ 



/■Spp/4 pi pi 

°H C (s P p) = dp 2 t dx p / dxp^2fi(x p , Pt)fj(xp, Pt) 

J pi min J m/ a pp J m/ x P s pp i,j 

d(Tij\XpXp-Spp) Pt) 



(16) 



dpt 

It is assumed, (as in ref.[[0],@]]) that secondary hard scatters are independent 
of each other so that the probability for the number of hard scatters h, in a 
pp collision at a given value of impact parameter, is given by the Poissonian 
probability distribution, 

Ph = exp (-(n(b, s P p))). (17) 

The inelastic hard cross section for pp— >-partons with Pt>Ptm%n can be written 
as : 



/oo 
db 2 p h 
h=i 



= nj db 2 [l - exp(-(n(6, s P p)))]. (18) 

The hard multiparton part || of the model was already implemented as 
part of the HERWIG Monte Carlo program . Here we give a brief summary. 

Event simulation starts with HERWIG generating a hard subprocess ac- 
cording to the leading order cross section. Both incoming and outgoing 
partons involved in the hard subprocess (at some hard scale Q) are evolved 
(backward for the incoming parton) through a coherent parton shower al- 
gorithm until they reach a typical hadronic scale, Qh w lGeV. The coher- 
ent parton shower algorithm [TIJ resums to all orders both single logarith- 



mic terms associated with the collinear emission (through implementation 
of the DGLAP splitting function) and single and double logarithmic terms 
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associated with soft emissions. All the successive partons emitted in the 
parton shower are colour connected in such a way that partons carrying 
colour (antico lour) tend to end up close in momentum and real-space to their 
anticolour(colour) partners. Once the evolved partons have reached a hadro- 
nisation scale Qh, clusters of partons are formed where each colour connected 
pair of partons forms a single colourless cluster. Each cluster is then de- 



cayed into hadrons according to phase space arguments ||11|| . Once the main 
hard subprocess has been simulated, energy-momentum conservation is used 
to calculate the momentum of the remaining diquark(antidiquark) or pro- 
ton(antiproton) remnant. Each diquark carries opposite colour to that of 
the parton involved in the hard subrocess. After the first hard scattering has 
been completed a number of secondary hard subprocess are simulated where 
each coloured remnant from the previous interaction is labeled as the new 
incoming hadron from which the quark valence distribution functions have 
been removed. 

Comparing the eikonal expression in (10) with the cross section in (18), 
we define the hard part of the eikonal as : 

XQCD(b, s) = ~(n(b, s P p)). (19) 

The QCD perturbative 2-to-2 parton-parton differential cross section di- 
verges as the transverse momentum of the scattering, p t , goes to zero. One 
must then fix a minimum value Ptmin for pt large enough so that the resulting 
cross section is not larger than the total inelastic non-diffractive cross section 
and also large enough for perturbative QCD to be reliable. This dependence 
of the perturbative QCD differential cross section implies that the eikonalised 
inelastic hard cross section would also depend on the arbitrary value chosen 
for the minimum p t . 



3.2 The expression for Xtotai{b,s) 

Our goal is to formulate a model which will be to a certain degree independent 
of Ptmin cutoff. To this end, to XQCD(b, s P p) which describes hard interactions 
with pt> Ptmin, we will add Xsoft(b, s P p), which will describe interactions with 
< pt<Ptmin- The full eikonal can now be expressed as : 



Xtotalip, S) = XQCD^b, S pf ) + Xsoftip, Spp). (20) 

To specify Xsoft(b, s P p) we will use a model introduced by Chou and Yang 
I2"f and Durand and Pi ||, which postulates that the elastic scattering is 



the shadow of the absorption resulting from the passage of one hadronic mass 
distribution through another. The transverse distribution of the matter is 
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assumed then to have the same shape as the charge distribution, as measured 
by the electromagnetic form factor, so that Xsoft{b, s) is of the form: 

Xsoft = C( Spp )A(b), (21) 

where C(s P p) is a constant to be determined below and A(b) is given in (13). 
The eikonal now consists of two parts, the hard, 

Xqcd = ~aT( Spf )A(b), (22) 

and the soft, 

Xsoft = ^softMAQ). (23) 

If we now consider the total cross section (8) with x(b, s)=Xtotai(b, s) and 
assume Xtotai(b,s) to be small, we can expand the exponential in (8). Inte- 
grating over the impact parameter b we get: 

&total(Spp) ~ &sOFt( S Pp) ®H i s pp) (^4) 

Following this approximation we will assume that our &soft( s pp) * s a Dare 
non-perturbative cross section for soft proton-antiproton interaction. Since 
&soft( s pp)> the bare soft cross section in (23), is not directly calculable in 
our model, its value can be determined from the experimental data, i.e. by 



using the total cross section measured by the CDF collaboration [13|. 



To give a prediction for future LHC data we will rely on the Donnachie- 



Landshoff model of the Pomeron [14] which fits successfully all the experi- 



mental data on total cross section for proton-proton scattering, 

/ s \ a p(°) _1 

CTtot(Spp) = crtot(so) — , (25) 

\s / 

where a p (t) = 1 + e + a p t with e = 0.08 and a' p = 0.25 GeV~ 2 . Assuming 
&tot{s PP ) = crtot(spp) at high energy, the total cross section at the LHC can 
then be written as 

/ , \ 2x0.08 



Vtot(Spp) = CTtot(scDF) yj^j ■ ( 26 ) 

We consider the value of the total cross section as the only phenomenolog- 
ical input to our model (in addition to those that are already in HERWIG0 
and the multiparton hard model M, which we leave at their default values). 

In addition the model contains the following assumptions : 
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• there can be more than one soft interaction in a soft proton-antiproton 
interaction. 



• the probability distribution for having m and only m soft scatters in 
a given proton-antiproton interaction P m obeys the Poissonian distribution 
(i.e. multiple soft scatters are uncorrelated) , 

Pm= {2Xs ° f . t)m exp(-2 Xsoft ). (27) 
ml 

• the hard and soft scatters are assumed to be independent i.e the prob- 
ability distribution for having h hard, and m soft, scatters is given by Phm- 



Ph, m {b, s P p) = n 1 ex P {-^XtotaiJ- (28 

hi ml 



We can now express cross sections as: 



a inel (s pf )=7r db 2 ]T P h , m (b,s pP ) = 7c ( ft2(l-e-WV l5 )) | ( 29 ) 

J h+m>l J ° 

the inelastic cross section with at least one hard scattering: 

<7 h ,inel(Spp) =71 [db 2 P h,Ub, Spp) = * T db 2 (l - e^W), 

J fc>l,m>0 

(30) 

and the total cross section: 

POO 

cT to t(s pf ) = 2tt ^ db 2 (l - e^W). (31) 

The cr^Q FT (spp) in (23) is then determined by fixing its value so that the 
the total cross section (31), at some ptmin, is equal to the total cross section 
measured by the CDF collaboration (a to tai cdf = 81.8±2.3mb) (see Sec. 5 for 
different values of &soft( s pp) an d Ptmin used). 

In order to have a realistic Monte Carlo event generator we have imple- 
mented our eikonal model as part of the HERWIG Monte Carlo program |3|] 
with the full use of its hadronisation model and parton showering properties. 

4 The Monte Carlo implementation of hard 
and soft processes 

In the following sections we will describe soft processes in detail. The hard 
processes are as in || (see Sec. 3). Before we start with the simulation of an 
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event we need to fix <Jsoft(, s pp)i the bare soft cross section (as explained in 
Sec. 3). Then by drawing from the probability distribution (28) we decide 
the number of soft and hard interactions in a given event. The event simu- 
lation can then start with HERWIG generating a hard subprocess followed 
by multiparton ones (see Sec. 3). Once the multiple hard scatters have been 
exhausted the soft multiple scatters are generated. 

4.1 Assumptions behind the soft subprocess 

At high s P p (center of mass energy) and low Q scale of interaction, i.e low-x, 
we might expect a proliferation of soft gluons. In order to simulate each soft 
subprocess we assume that the bare soft cross section & 1 soft( s pp)-i corresponds 
to a soft interaction between two incoming soft gluons. That is, we model 
the soft inelastic collision between two remnants as a soft elastic collision 
between two partons within them (see fig.l). 

All soft gluons carry a colour charge and have an effective mass m g (in- 
spired by HERWIG's model of the hadronisation of the outgoing gluons). 
The two outgoing effective gluons in fig.l are on mass shell and colour con- 
nected to the remnants and to each other, as shown in fig.l, each gluon 
colour(anticolour) is connected to its anticolour(colour) partner. 



Figure 1: &soft( s pp) corresponds to a soft collision between the two soft gluons (full 
color picture). Remnants are also connected to each other via t channel gluon line. 

Because we wanted to keep the implementation of our model as simple 
as possible we have decided to sever the colour connections between the 
remnants and the outgoing soft gluons (see fig. 2, dashed solid lines show 
severed colour connections). We expect this to be a reasonable approximation 
since the largest part of the phase space, in the gap between the two remnants, 
will be filled by final state hadrons produced by the colour field stretched in 
between the two outgoing soft gluons (note that in both figs 1 and 2, two 
colour lines stretch across the central region). In the hadronisation phase 



P-Remnl 



2 




J^remnl 



Pbar-remn2 
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Figure 2: The soft collision between the two soft gluons, with dashed solid lines indicating 
the severed color connections between the remnants and the outgoing gluons, forming two 
clusters, q\qi and qiqi. 

the two outgoing soft gluons form two clusters (qiq 2 , Q2Q1) (see fig. 2) which 
are decayed into final state hadrons (using HERWIG's cluster hadronisation 
model @). 

4.2 Implementation of the soft process 

Following fig. 2, for each soft remnant-remnant interaction we generate two 
soft effective gluons, on mass shell, with previously tuned masses m g = 0.75 
GeV used in the HERWIG hadronisation model || . For each soft interaction 
we will go through the following chain of events; first the maximum fractional 
longitudinal momentum allowed for each soft gluon (x g i umax i j2 ) is determined 
using the light-cone definition of the longitudinal momentum fraction x of 
the two remnants, 

_ _ ^remnl,2 "I" Pz r emni,2 (*iO\ 

•Eglumax 1,2 •^remn 1,2 j-, j p ■ v" / 

The longitudinal momentum fractions for each soft gluon is then sampled 
from an f (x) sea P artons = - distribution (which we consider reasonable since 
we expect the effective gluons to be Regge-like) between some minimum value 
Xmin (the cutoff) and the maximum value, x g i umax i^, allowed by (32). Once 
their fractional longitudinal momentum has been determined, a transverse 
momentum for each soft gluon is then sampled from a Gaussian distribution 
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for values of p t , 0<Pt<Ptmin, according to 



dNs °!^ luons = D exp < p t < p tmin . (33) 

Before sampling the transverse momentum of each soft gluon, we will have to 
determine the slope (3 of the particle p t distribution in the central region and 
the normalisation constant D in (33). We impose then another condition on 
the p t distribution, namely, that the transverse momentum distribution in p t 
of soft and hard gluons should be continuous at the p t cutoff (34) (the same 
procedure was implemented in fll5l). 

A AT. 

(34) 



dNsoft gluons i _ G^-^/iard gluons i 

C[2p t \pt=Ptmin ~ d^Pt VPt=Ptmin- 



We have then to solve two equations with two unknowns. The first condition 
is that the number of soft gluons should correspond to the soft cross section 

a SOFT( s Pp) : 

[ Ptmin d 2 p t D exp {-(3p 2 t ) = 4oV(%0, (35) 

J 

the second condition is that of the smooth transition between the soft and 
hard perturbative gluons' transverse momentum distributions at the value of 
p t cutoff (34), re-expressed as: 

r\ I T3 2 \ d(j(s p p, ptmin) foc\ 

D exp {-f3p tmin ) = t~2 > ( 36 ) 



dp, 



tmi'n 



w j iere da ~( s p P ^Ptrmn) j g f u }} differential parton-parton hard cross section cal- 

" train 

culated at the value of p t = Ptmin- 

In order to determine the four momentum of each outgoing remnant and 
each soft gluon, two additional conditions need to be satisfied simultaneously; 
that of the total energy-momentum conservation between the initial beam- 
beam remnants (remn 12 ) and final (remn 12 + soft gluons) states and that 
the outgoing soft gluons and remnants are on mass shell, with each remnant 
having the same mass before and after the soft interaction. 

Once all four momenta of the outgoing remnants and soft gluons have 
been determined a new maximum longitudinal fractional momentum x g i umax \^ 
is calculated, according to (32), for the next pair of soft gluons. The same 
chain of events, as described above, is then iterated until all soft interac- 
tions have been generated. If energy-momentum conservation is violated the 
scattering is vetoed. 
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5 Comparison between the eikonal model and 
data from the CDF collaboration 



We will compare our model predictions with the experimental data coming 
from the study of the behaviour of the underlying event in hard scatter- 
ing proton- antiproton collisions at 1.8 TeV. In ref. |I[ the underlying event 
in hard proton-antiproton collisions is defined as all the ambient interac- 
tions surrounding the hard scattered jets. These ambient interactions include 
beam-beam remnants, initial and final-state radiation and multiple hard in- 
teractions, if they do occur. The beam-beam remnants are just what is left 
over after a parton has been taken out from each of the initial incoming 
hadrons. The jets in a hard interaction are produced by streams of outgoing 
hadrons with large transverse momentum, from partons that have undergone 
a hard two-to-two scattering. In interactions where a hard two-to-two sub- 
process occurred, jets of hadrons are produced approximately back to back in 
azimuthal angle. This simple jet structure allows for the study of the particle 
distribution in azimuthal angle <fi. There are three regions of study separated 
in A0 from the leading highest p t jet (see fig. 3). The "toward" region contain- 
ing the leading charged particle jets, the "away" region containing particles 
produced opposite in from the leading jet and the "transverse" region per- 
pendicular to the plane of the hard two-to-two scattering which is the most 
sensitive region to the underlying event. 



Jet#1 Direction 

A 




Figure 3: Toward, away and transverse regions from the leading inp t jet (see ref. [Q).The 
angle A(f> — 4> - 4>j e t#x is the relative azimuthal angle between charged particles and the 
direction of jet#l. The three regions are as follows; the toward region defined over the 
range |A</>| < 60° (this region includes the particles of the leading jet#l), the away region 
defined over the range |A(/>| > 120° and the transverse region defined over the range 
60° < |A0| < 120°. Each region, toward, transverse, and away covers the same range 
|Ar ? |x|A(/.|=2xl20 o . 
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The simple jet algorithm used in builds jets from charged particles with 
p t > 0.5 GeV and \r]\ < 1 and is restricted to particle jets with transverse 
momentum less than 50 GeV. The charged particles form jets of radius R 
= 0.7 in r)-(j) space which contain particles from both underlying event and 
hard scattering. Every charged particle is assigned to a jet with a possibility 
that some jets might consist of just one charged particle. The transverse 
momentum of the jet is defined as the scalar p t sum of all particles within 
the jet. Once all jets have been produced the one with the highest p t (jet#l) 
is defined to be the leading jet. The experimental data are uncorrected and 
the theoretical simulated data are corrected for the track finding efficiency 
by removing, on average 8%, of the charged particles. 

5.1 Results 

A comparison between simulated (corrected) and experimental data is made 
and presented in figs 4-7, with ptmin = 2.5 GeV. The agreement seen in 
figs. 6 and 7 is particularly important since the transverse region is the most 
sensitive to the underlying event activity. In figs. 8 and 9 we compare our 
predictions for the transverse region with those of the existing models (with- 
out including Pythia), HERWIG + Underlying Event model @, HERWIG + 
Multiparton Hard model In HERWIG's default Underlying Event model, 
the activity is enhanced in the transverse region by adding to its hard part 
a Minimum Bias Event based on the parametrisation of the UA5 pp 
minimum-bias data. As shown in figs. 8 and 9 and in |l| this model fails to 
produce enough activity for both the average number of charged particles 
and the average scalar p t sum and in addition has a wrong (too steep) p t 
dependence. Overall our eikonal model is in better agreement with the ex- 
perimental data than either of the two underlying event models. This result 
is even more significant since we did not fit to the experimental data but 
rather, make predictions based on the value of the total cross section used 
as phenomenological input to our model. 

5.2 The invar iance of the model to ptmin 

In the multiparton scattering model Q , the parameter p tm i n is a cutoff scale 
and plays a crucial role in determining the prediction of the model. In our 
eikonal model it is rather a matching scale between the hard and soft parts 
of the model. If the matching works perfectly, the results should be ptmin 
independent. 

To test the invariance of the eikonal model to its Ptmin parameter we 
simulate two sets of data, with two different values of Ptmin = 2 GeV and 3 



12 



GeV in addition to the value of 2.5 GeV, which we have already shown. The 
values of inclusive cross sections with the average number of soft and hard 
scatters are presented in table. 1. 



Ptmin (GeV) 


Vsoft(spp) (mb) 


ar( Sp p)(mb) 


(n S oft) 


( n hard) 


2.0 


39.7 


99.2 


0.7 


1.7 


2.5 


85.6 


51.3 


1.5 


0.9 


3.0 


109.7 


28.7 


1.9 


0.5 



Table 1: Values of inclusive cross sections and average numbers of hard and 
soft scatters. 

We show the results in fig. 10. We see that as expected the results in 
the multiparton model are strongly ptmin dependent, with smaller values of 
Ptmin producing more activity. However, even the smallest value of ptmin is 
well below the data and no value of ptmin for which the hard part of HER- 
WIG is reliable gives a good description of the data (although it is worth 
noting that with ptmin fixed at 2.5 GeV it is possible to fit the data by in- 
creasing the proton radius by a factor of about 1.7 fL7| . Since this model 
does not provide any explanation of how the radius should vary as a function 
of center-of-mass energy, or for different collision types, we do not consider 
it to be useful phenomenologically) . In contrast, our eikonal model has a 
smaller ptmin dependence, with the additional soft scatters providing extra 
activity to compensate for that lost from the hard scatters as ptmin increases. 
However, we see that they actually overcompensate, with decreasing ptmin 
leading to decreasing activity. This residual ptmin dependence indicates that 
the matching of the soft and the hard parts of the model is still not perfect 
and can perhaps be improved with further refinements of the model. Nev- 
ertheless, we are satisfied that with such a simple, physically motivated and 
parameter free model, we have provided a significant improvement to the 
Ptmin dependence and the description of data. 



13 




10 15 20 25 30 35 40 45 50 

Pt leading jet(GeV) 



Fi gure 4: The average Pt sum of charged particles as a function of Pt (leading charged jet) 
in the toward region. HERWIG + Eikonal model solid line simulated data, experimental 
data solid circles. 
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Figure 5: The average Pt sum of charged particles as a function of Pt (leading charged 
jet) in the away region. HERWIG + Eikonal model solid line simulated data, experimental 
data solid circles. 
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Figure 6: The average Pt sum of charged particles as a function of Pt (leading charged 
jet) in the transverse region. HERWIG + Eikonal model solid line simulated data, exper- 
imental data solid circles. 
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Fi guxe 7; The average number of charged particles as a function of Pt (leading charged 
jet) in the transverse region. HERWIG + Eikonal model solid line simulated data, exper- 
imental data solid circles. 
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Figure 8: The average Pt sum of charged particles as a function of Pt (leading charged jet) 
in the transverse region. HERWIG + Eikonal model (solid line), HERWIG + Underlying 
Event model (solid dashed), HERWIG + Multiparton Hard model (dotted), experimental 
data solid circles. 
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Figure 9: The average number of charged particles as a function of Pt (leading charged 
jet) in the transverse region. HERWIG + Eikonal model (solid line), HERWIG + Under- 
lying Event model (solid dashed), HERWIG + Multiparton Hard model (dotted), experi- 
mental data solid circles. 
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10 15 20 25 30 35 40 45 50 

Pt leading jet(GeV) 



Figure 10: The average Pt sum of charged particles as a function of Pt (leading charged 
jet) in the transverse region. HERWIG + Eikonal Model for the two sets of ptmm=3.0 
GeV (solid line) and 2.0 GeV (dashed), HERWIG + Multiparton Hard Model p tmin =2.0 
GeV (dotted), p tm i„=3.0 GeV (dott dashed), experimental data points. 
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6 Discussion and Conclusion 



In this paper we have proposed a Monte Carlo model, based on a simple 
eikonal model, for multiparticle production in hard proton-antiproton col- 
lisions. The model contains two parts, the hard part, which was already 
implemented [|j] (running in conjunction with HERWIG), to which we have 
added a soft part, which uses showering and hadronisation models from HER- 
WIG and allows us to extend our simulations to the non-perturbative soft 
region (i.e. particles with < p t < Ptmin)- One of our goals was to pro- 
duce a simple model (with a minimum number of parameters) which would 
be reasonably independent of the value chosen for the ptmin- Furthermore, 
the eikonal multiparticle approach to hadron-hardon scattering is particu- 
larly interesting since it prevents the calculated cross sections from violating 
unitarity at high energies. By fixing one parameter in our model, namely 
the total cross section, we provide a good description of the measured data 
from with a reasonably small ptmin dependence. We have also shown that 
our model gives a better prediction than either the HERWIG + Underlying 
Event or HERWIG + Hard Multiparton @] models. 

Our ultimate goal is to determine how much of the underlying event activ- 
ity in the transverse region is due to the perturbative physics (by including 
NLO calculations for example) and how much is due to non-perturbative 
physics, and to make more reliable predictions of the underlying event. We 
are planning to test the model with the lower energy UAl collaboration data, 
HERA data, higher p t CDF data and give a prediction for the underlying 
event activity at the LHC. 
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